Quantum fluctuations of a vortex in an optical lattice 
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Using a variational ansatz for the wave function of the Bose-Einstein condensate, we develop a 
quantum theory of vortices and quadrupole modes in a one-dimensional optical lattice. We study the 
coupling between the quadrupole modes and Kelvin modes, which turns out to be formally analogous 
to the theory of parametric processes in quantum optics. This leads to the possibility of squeezing 
vortices. We solve the quantum multimode problem for the Kelvin modes and quadrupole modes 
numerically and find properties that cannot be explained with a simple linear-response theory. 
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Introduction. — There is a large body of research on 
vortices in Bose-Einstein condensates, both before and 
after their first experimental observation Q, Q. A re- 
cent revievif of this vibrant field was given by Fetter and 
Svidzinsky One reason for this interest is that vor- 
tices play a crucial role in explaining rotational and dissi- 
pative properties of superfluids. They are also important 
for an understanding of classical fluids and wave excita- 
tions of a vortex were therefore already studied by Lord 
Kelvin. Such excitations have been observed in super- 
fluid helium 4\ and, most recently, in Bose-Einstein con- 
densates by Bretin et al. j^. The latter experiment also 
elucidated the role of the quadrupole modes of a Bose- 
Einstein condensate in exciting the Kelvin modes of the 
vortex line. A theoretical understanding of the observed 
phenomena was provided by Mizushima et al. 0. 

The phenomena of superfluidity is closely connected 
to the existence of quantized vortices and relies on phase 
coherence across the condensate. Placing the condensate 
in an optical lattice makes the establishment of phase co- 
herence more difficult and can even induce a transition 
into a Mott-insulator state. This quantum phase tran- 
sition was indeed recently observed experimentally 
In a one-dimensional optical lattice the condensate splits 
into a stack of weakly-coupled disk-shaped condensates, 
which leads to some intriguing analogs with high- Tc su- 
perconductors due to their similar layered structure 

In this Letter we combine the above ideas and 
study a vortex in a Bose-Einstein condensate in a one- 
dimensional optical lattice. This problem has not yet 
been addressed either experimentally or theoretically, but 
experimental studies are possible with current techni- 
cal capabilities. Indeed such an experimental setup ap- 
pears to be a promising way to achieve the quantum Hall 
regime in a Bose gas, for which we need to have about 
one vortex per particle 01 ■ Also outside the quantum Hall 
regime an optical lattice allows for unique possibilities to 
study vortices and their quantum fluctuations. In a nor- 
mal three-dimensional condensate a vortex line is a very 
complicated dynamical object and it is nearly impossible 
to tune its properties in a controlled way. Placing the 
vortex in an optical lattice provides a way, for example. 



to reduce the stiffness of the the vortex line by increasing 
the depth of the lattice. 

With this goal in mind, we develop a quantum the- 
ory for the Kelvin modes in an optical lattice and ob- 
tain the Kelvin mode dispersion relation analytically by 
means of a variational approach. Furthermore, because 
of their importance for the creation of the Kelvin waves, 
we study in detail the coupling between the Kelvin modes 
and quadrupole modes. Our approach leads to a quan- 
tum multimode problem with an interaction between 
Kelvin modes and quadrupole modes that closely resem- 
bles squeezing Hamiltonians familiar from quantum op- 
tics. We solve this quantum multimode problem and 
show how the squeezing is reflected in the quantum me- 
chanical uncertainty of the vortex line. In addition, nu- 
merical solution of the multimode problem reveals prop- 
erties of the nonlinear dynamics that cannot be explained 
with a linear-response theory. 

Theory. — The theory we use is based in our previous 
work |lfllll| . and we briefly summarize the main points. 
We consider an elongated Bose-Einstein condensate in a 
magnetic trap with a radial trapping frequency tOr that 
is much larger than the axial trapping frequency w^. In 
addition to this cigar-shaped potential the condensate 
experiences an one-dimensional optical lattice Vq (r) = 
Vo sin^ (27rz/A), where Vq is the lattice depth and A is the 
wavelength of the laser light. Furthermore, we assume 
that the lattice is deep enough so that it dominates over 
the magnetic trapping potential in the z direction and 
that the number of sites are very large. Then we can in 
first instance ignore the magnetic trap in the z direction. 

We expand the condensate wave function in terms of 
wave functions (a;, y) $(z) that are well localized in 
the sites labeled by n. We do not yet assume any- 
thing about the wave functions <&„ {x, y) for the two- 
dimensional condensates, but for the wave function in the 
z direction, $ (z), we use the ground-state wave function 
of the harmonic approximation to the lattice potential 
near the lattice minimum. This harmonic trap has the 
frequency ujl — (27r/A) A/2Vb/m, where m is the atomic 
mass. These considerations lead to an energy functional 
of the weakly-coupled two-dimensional condensates. 
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which is characterized by the strength of the on-site par- 
ticle interactions Anil = AttN (a/lr) \/ uoj^jl-KUjir and by 
the strength of the nearest-neighbor Josephson coupling 
J = (I/Stt^) (Xjl^f [7rV4 - l] e-(V4;.)^ ^ . 

In these formulae N is the number of atoms in each site, 
a is the scattering length, and Z^, = ^JhJmujL. Also, we 
use trap units so that the unit of energy is hujri the unit 
of time is l/wr, and the unit of length is Ir = ^hjmLOr- 
Our theory is valid when the lattice depth is much big- 
ger than the chemical potential. Also, the strength of 
the Josephson coupling should be small, but still large 
enough to support a superfluid state. The superfluid 
Mott-insulator transition in one-dimensional optical lat- 
tice is expected to occur only for very deep lattices |0, 
and such a transition has indeed not yet been observed 



experimentally. 

Variational ansatz. — A variational ansatz for the two- 
dimensional condensates should satisfy two obvious phys- 
ical criteria. First, the ansatz must predict the vortex 
precession dynamics correctly. Second, the ansatz must 
predict the quadrupole mode frequencies with sufficient 
accuracy. In an earlier publication we used an ansatz 
that fulfilled the second requirement very well in the full 
parameter range from the non-interacting limit to the 
Thomas- Fermi regime Unfortunately, this ansatz 

does not fulfill the first requirement because it does not 
properly account for the size of the vortex core, which is 
not important for the collective modes but plays a crucial 
role in the dynamics of the vortex. Therefore, we use the 
ansatz 



{x,y) cx exp ( 



(x^ - y'^) - exv,n{t)xy + i tan 



X - Xnit) 



(1) 



where (xn, y-n) is the position of the vortex in the nth site 
and the variational parameters e„ and txy,n are complex. 
In addition, the square of the wave function is normalized 
to N . We have chosen our variational parameters such 
that only Kelvin modes and quadrupole modes are in- 
cluded. Our theory can be easily adapted to include the 
monopole mode as well, but we choose to focus on the 
quadrupole modes due to their experimental importance 
for the creation of Kelvin waves 0. Since our ansatz 
does not have a core, the average kinetic energy diverges. 
Therefore, we introduce a small distance cut-off ^ as an 
additional variational parameter. Together with this 
parameter will be determined by minimizing the equilib- 
rium energy functional. This minimization is a simple 
numerical task that can be done once in the beginning 
of the calculation. In that manner the condensate with 
a vortex is slightly bigger than a condensate without a 
vortex and the vortex core size becomes smaller with in- 
creasing atom numbers. Quantitatively, we find to a very 
good accuracy that ^ = a/Bq. 

The ansatz in Eq. ^ fulfills both our requirements 
with physically realistic parameter values and is still sim- 
ple enough to be analytically tractable. Our ansatz fails 
for very weakly-interacting condensates, but provides an 
accurate description when the interaction strength U is 
larger than about 10. In the other regime our previous 
ansatz can be used |0 . In this Letter we do not consider 
rotation of the trap, but a change to the rotating frame 
can again be easily included into our approach. 

Quantization of the vortex and the quadrupole 
modes. — If we expand the Lagrangian for our system 
to second order in the deviations from equilibrium, we 



are in a position to study the eigenmodes of the system. 
The eigenmodes we find are the Kelvin modes and two 
quadrupole modes with quantum numbers m = ±2. To 
second order there is no coupling between the vortex po- 
sitions and the quadrupole modes and we can treat the 
two independently. We transform to momentum space 
using the convention fk — l/v^/Vg e^*'^^"/,!, where 
Ns is the number of sites. We then proceed by defin- 
ing the vortex positions in terms of bosonic creation and 
annihilation operators as Xk = {vlk+Vk)/2VNB^ and 
Vk = i{vlk - Vk)/2y/NB^ [13. With this definition the 
Lagrangian for the vortex position takes the desired form 

Lk^Y^ (ivivk - i^K (fc) vivk^ , (2) 
fc 

where, to leading order in Bq, the Kelvin mode dispersion 
is given by (fc) = Bo/2 + T [O, B^] (2J{k) - Bo/2). In 
this expression J{k) = J (1 — cos (/cd)), d = A/2, and 
r [a, z] is the incomplete gamma function. Our result for 
the Kelvin mode dispersion relation has some interesting 
features as seen in Fig. ^ (b). First, the initially nega- 
tive precession frequency, that indicates that a straight 
vortex line is unstable without rotation, can change sign 
if the momentum k of the Kelvin mode is large enough. 
Moreover, the dispersion is quadratic for small momenta. 
Both features are due to the attractive interactions be- 
tween neighboring pancake vortices, which is harmonic 
for small separations and due to the energy cost for 
having phase differences between sites. The quadratic 
Kelvin dispersion in a lattice contrasts with the behavior 
of the vortex line in a three-dimensional bulk superfluid. 
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where the Kelvin mode dispersion relation behaves like 
fc^ln(l/A;^) and has a logarithmic dependence for small 
k. Third, the vortex position is "smeared" by quantum 
fluctuations. In particular we have for the Kelvin mode 
vacuum state that (x^) = (yl) = l/(47VBo). There- 
fore, the quantum properties of the vortex become more 
important in a lattice. This is due to the reduced on- 



site particle number, as opposed to the total number of 
particles, and the spreading out of the condensate wave 
function as the lattice depth is increased. 

The quadrupole modes are technically somewhat more 
complicated, but can be quantized in a similar way as 
the Kelvin modes. The quadrupole mode frequencies are 



J 



UBo 



SoT [0,^2] 



2J{k) [A{k) 



A{k) 



±Bo, 



(3) 



r 



where 



Aik) 



8BoJ{k) 



2T [0,^2] 



^ 4 + 8Bo./(fc) + (1 
The Lagrangian for the quadrupole modes is then 



EE 

Ti=±2 k 



kQm,k - UJm (fc) qJn k^m.k 



(4) 



(5) 



In Fig. n (a) we plot the behavior of the mode frequen- 
cies at A: = as a function of the interaction strength. 
The quadrupole mode frequencies are in good agreement 
with the ones calculated using the Bogoliubov-de Gennes 
equation and with realistic interaction strengths the 
vortex precession frequency is close to the result using 
a Thomas-Fcrmi wave-function for the disk-shaped con- 
densate |3|. 

(a) (b) 
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FIG. 1: (a) Frequencies of the m — ±2 quadrupole modes 
and the Kelvin mode mode at = as a function of inter- 
action strength, (b) The constant line is LO-2 (0) /2 and the 
curve is the Kelvin mode dispersion. In linear response the 
energy transfer between the m — —2 quadrupole modes and 
the Kelvin modes becomes nonvanishing if the Kelvin mode 
frequency matches lu-2 (0) /2, i.e., when the lines cross. We 
used U = 100 and J = 0.1. 

Coupled Kelvin modes and quadrupole modes. — We 
now proceed to apply our previous results to study the 
coupling between the Kelvin modes and the quadrupole 
modes. To obtain coupling between the quadrupole 
modes and the Kelvin modes, we must expand the La- 
grangian up to third order. At third order there are con- 
tributions from the kinetic energy, the time-derivative 



term in the Lagrangian, and the Josephson coupling 
term. Since we assume small J we can safely ignore 
the last one. The remaining terms describe the coupling 
between the Kelvin modes and both quadrupole modes. 
The contribution from the time-derivative term in the 
Lagrangian includes the time-derivative of the Kelvin 
mode operators. In the rotating-wave approximation we 
can replace these time-derivatives with the Kelvin mode 
dispersion. In addition, this eliminates the nonresonant 
terms from the Hamiltonian, which in particular removes 
the coupling between the Kelvin modes and the m = 2 
quadrupole mode. In this way we obtain the following 
interaction Hamiltonian for the Kelvin modes and the 
quadrupole modes 

Hi ^'^Ec (k, k') {q-2,k'^k+k'^^-k + q^-^2.k'^-k-k'Vk 

k,k' 

(6) 

where E,{k,k') = ^A{k')/8Ns (BqF [O, B^] +ujK{k)). 
This expression is accurate as long as the coupling be- 
tween the Kelvin modes and the quadrupole modes is 
sharply peaked around some resonant value kg. Note 
that the angular momentum is conserved, since the 
quadrupole mode has an angular momentum of —2 
whereas each Kelvin mode has an angular momentum 
of-1. 

Multimode problem. — Our interaction Hamiltonian 
shows how the coupling singles out a wavenumber for 
the Kelvin modes and thus selects the wavelength for the 
ensuing wiggles in the vortex line We demonstrate 
this in Fig. ^ (b). The width of the resonance is still 
unknown, but we can give a simple estimate. Assume 
that initially we have a straight vortex line and a coher- 
ent quadrupole field (g_2,fc(0)) — g-2,o(0)(5fe,o- At short 
times we can assume that the quadrupole field stays con- 
stant and then solve the operator equation for the Kelvin 
modes v±k. It turns out that the solutions are dynami- 
cally unstable only when |Aa;(A;) | < |£^c (fc, 0) g_2,o(0)|, 
where AcL'(fc) = uj-2 (0) — 2uik {k). Our numerical solu- 
tions of the multimode problem confirm this simple esti- 
mate with a reasonable accuracy. The number of modes 
within the interval depends on the number of sites in the 
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lattice and on the quadrupole amplitude, but with real- 
istic parameter values the multimode behavior is more 
likely than a single-mode behavior. 

In linear response we can use Fermi's golden rule to 
calculate the transition rate from the initially coherent 
state IV'i) = \q-2,o,{'^k}) = \q-2,o,{^})- After a short 
calculation we get the transition rate as 

^ 7rE,{ko,0)^\q.2,o (0)P . . 

'^^^ 2jr[0,B2]|sin(fcod)| ' ^ ^ 

In Fig.|2]we show an example of the time evolution of the 
system, when the quadrupole field is treated as a clas- 
sical field, but Kelvin modes are treated fully quantum 
mechanically. In the same figure, we also compare the ex- 
ponential g_2,o (0) e"'"''^^* with the numerical solution 
of the nonlinear multimode problem and show the num- 
ber distribution of Kelvin modes at long times. It should 
be noted that at these time scales the system appears 
irreversible due to the large number of modes present. 
If the number of sites in the lattice is smaller the os- 
cillatory behavior becomes stronger. It should also be 
realized that, depending on parameters, the result from 
Fermi's golden rule can be unreliable since it ignores the 
finite lifetime of the final states. It is interesting to ob- 
serve that the number distribution of the Kelvin modes 
is asymmetric. This is due to the asymmetric nature of 
the coupling constant Ec (fc, 0) and the detuning Auj (fc). 

The interaction Hamiltonian in Eq. © closely resem- 
bles the squeezing Hamiltonian encountered for paramet- 
ric processes in quantum-optics [l^ . Therefore, it is not 
surprising that we can observe a dramatic squeezing un- 
der certain conditions. For example, if we only have a sin- 
gle resonant mode then the initial minimum uncertainty 
state, i.e., the vacuum, will remain a minimum uncer- 
tainty state, but can become strongly squeezed. Such a 
squeezing is reflected in the vortex position distribution. 
This is shown in Fig. |21 where we also monitor the time- 
evolution of the deformation e(t) — max ;^?l!l!T/^!i 
of the vortex position distribution. Because the uncer- 
tainty ellipse that characterizes the vortex position dis- 
tribution precesses in time, the deformation is maximized 
with respect to the rotation angle of the coordinate sys- 
tem. 

Discussion. — We have developed a quantum theory 
of a vortex in a one-dimensional optical lattice and ap- 
plied the theory to study the coupled quadrupole-Kelvin 
mode dynamics. Our solution of the Kelvin mode dis- 
persion relation essentially amounts to a solution of the 
vortex-line dynamics in an optical lattice. Among other 
things our study revealed the possibility of squeezed vor- 
tex states and asymmetric number distributions for the 
Kelvin modes. As yet, there have been no experiments 
on vortices in optical lattices reported. Detailed studies 
of the Kelvin mode dynamics would require imaging of 
the three-dimensional vortex line which is difficult, but 



not impossible |5|,ll5j. An alternative, but less direct way 
to probe the system, would be to track the evolution of 
the quadrupole oscillations. If the resonance is narrow 
enough our theory predicts revivals of these oscillations. 
In addition, the multimode nature of the Kelvin mode 
distribution can cause a decay rate of the quadrupole 
mode that is quite different from Fermi's golden rule. 
Our theory can be applied to a variety of phenomenon. 
For example, the interplay between the superfluid Mott- 
insulator transition, the quantum melting of a vortex lat- 
tice, and the quantum Hall regime in an optical lattice is 
a very interesting topic for further research. 



(a) (b) (c) 




FIG. 2: (a) The time-evolution of the quadrupole mode (solid 
line) compared with the result based on Fermi's golden rule 
(dashed line), (b) The number distribution of the Kelvin 
modes at the end of the simulation, when the quadrupole 
amplitude has essentially vanished. The circles indicate the 
individual Kelvin modes included in the simulation and the 
solid line is an interpolation between the data points, (c) 
The time-evolution of the asymmetry factor e{t). We used 
U = 100, J = 0.1, Ns = 1001, N = 100, and (g-2,o(0)) = 25. 
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